Upper Bound on the Energy of Particles and Their Secondary Neutrinos 
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We derive an upper limit to the energy of nuclei accelerated via the Fermi mechanism in any 
relativistic Shockwave, driven by any astrophysical engine. This bound is accessible to current and 
upcoming ultra-high energy neutrino experiments. Detection of a single neutrino with energy above 
the upper limit would exclude all sites of shock acceleration, and imply physics beyond the Standard 
Model. We comment on the possibility that relativistic flows launched by supermassive black hole 
mergers are the source of the observed ultra-high energy cosmic rays. 

PACS numbers: 95.85.Ry, 98.70.Rz, 98.70.-f 



The origin of ultra-high energy cosmic rays (> 
10 185 eV; UHECRs) is unknown. Broadly speaking, they 
can be produced by shock acceleration via the Fermi 
mechanism [D-Q], direct acceleration by the electromag- 
netic fields of compact objects and the decays of hy- 
pothetical super-massive particles beyond the Standard 
Model ( "top-down" scenarios) 0, ■ Fermi acceleration 
is thought to operate in supernova remnants, producing 
the low energy CRs (< 10 15 5 eV) in the Galaxy. Sim- 
ilarly, the relativistic shocks of gamma-ray bursts and 
the flares of active galactic nuclei are leading prospective 
astrophysical engines for UHECR acceleration 0, Q . 

Whatever their source, UHECR acceleration must oc- 
cur in the local universe. If UHECRs are protons, they 
are attenuated on a scale of < 100 Mpc for energies 
> 5 x 10 19 eV by the Greisen-Zatsepin-Kuz'min (GZK; 
[ill ]) mechanism, producing secondary particles that de- 
cay to e , 7-rays, and UHE neutrinos. If UHECRs are 
nuclei, photodisintegration on the cosmic infrared back- 
ground suppresses the flux at a similar energy 
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Although a decrease in the UHECR flux above the 
GZK threshold has now been observed 13, Il4| . it is 
unknown whether this is evidence of UHECR attenua- 
tion, or instead an intrinsic limitation to the accelera- 
tion process in astrophysical engines. If the former, then 
secondary UHE neutrinos with energies above the GZK 
cutoff should propagate to us unimpeded from across the 
universe [HI, [3] • Since experimental searches for UHE 
neutrinos are sensitive up to ~ 10 26 eV [13], it is of con- 
siderable interest to understand if physical bounds can 
be placed on the acceleration process that limit the max- 
imum energy of the primary UHECR spectrum. 

Such a bound has been constructed for the case of di- 
rect electromagnetic acceleration p|, but no such bound 
has been formulated for the case of relativistic shock ac- 
celeration. Here, we show that the physics of the Fermi 
mechanism and synchrotron cooling, together with an ab- 
solute upper bound on the bolometric luminosity of any 
astrophysical engine, combine to produce an experimen- 
tally testable upper bound on the energy of the UHECR 
population in the universe. 



The Maximum Energy of Particles. — We consider par- 
ticle acceleration in a flow of total luminosity L and bulk 
Lorentz factor T. As in Ref. Q, we assume that the 
central engine is transient and stationary in the observer 
frame, and that the region of particle acceleration is a 
shell of radial extent Sr = 6r' /V = r/T 2 , and perpendic- 
ular extent r'^ = rj_, where r is the distance from the 
shell to the engine, and where primes denote quantities 
in the flow rest frame. The connection between Sr, r, 
and r follows from the fact that an unsteady flow of thin 
shell geometry traveling at c — 5v, with Sv <C c, broadens 
kincmatically to a thickness Sr re rSv/c. 

The comoving magnetic energy density in the shell is 
u' B = eB-k/(47r\I'r 2 r 2 /3c), where /? = v/c is the shock ve- 
locity, ee(< 1) measures the efficiency of magnetic field 
generation, and W(< 1) = Q/Att measures the beaming of 
the flow into solid angle Q. A necessary requirement for 



acceleration is that 0, 0] 2irfR' L 



< mm 



[Sr', 



where 



R' L = E'/{ZeB') is the Larmor radius, E' = E/T = 
j'mc 2 is the particle energy, Ze is the charge, B' is the 
magnetic field strength, and / is a constant of order unity. 



We first focus on the case with Sr' < 



and return to 



the physics of more highly beamed flows below. Assum- 
ing P re 1, which we confirm below, and combining these 
expressions, one obtains an upper bound on the acceler- 
ated particle energy 



< 
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However, there are two crucial limitations to Eq. ([1}. 
First, the same magnetic field that confines the particles 
to the accelerating region also causes radiative losses via 
synchrotron. Acceleration to E requires that the accel- 
eration time t' acc ~ R'l/c is less than the synchrotron 
cooling time, t' nch , which sets a lower limit on T Q: 



r > 



E 



synch 
2/5 



^f(Ze) 3 
9m 2 c 5 St 



1/5 



2e B L 



1/10 



(2) 



where m is the particle mass, St ~ r/(r 2 c) is the ob- 
served variability timescale of the system, and we have 
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assumed that v'/c = 1, where v' is the particle velocity 
in the rest frame of the flow. Because f^cc ^ R'li * ne con " 
dition t' acc < t' aynch implies that radiation reaction can be 
ignored in deriving Eq. ([2]) fljj ]. 

The second fundamental constraint on Eq. ([1]) follows 
from the fact that there is a maximum L attainable by 
any astrophysical engine [19 ]: 



5 /G ~ 3.63 x 10 59 ergs s" 



(3) 



equivalent to radiating the entire rest mass of a body in 
a single light travel time across its Schwarzschild radius. 
Lmux is only approached during the merging of binary 
BHs in gravity waves. We return below to whether L max 
in gravity waves can couple to electromagnetic fields, thus 
producing u' B oc e^L with e B > 0. Here, it is sufficient 
to note that L < L max for any astrophysical engine. 

The limits on T and L combine with Eq. ([IJ to produce 
a unique E max at a critical r cr i t . We find that 
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is the maximum energy of particles accelerated in rela- 
tivistic shocks driven by maximally luminous astrophysi- 
cal engines. To evaluate £ mM , note that the observed 
variability timescale St = t/T 2 , where t is the light- 
travel time between the engine and the acceleration re- 
gion. Since t is bounded by the time for GZK losses for 
protons (fioss), or photodisintegration losses for nuclei, 



t (f)«lxl0 24 eV^ 4eB 
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and T crit « 250Z(e B // 2 / 3 *) 3 / 10 (™ p /m) 4 / 5 (fi O ssA) 1 / 5 , 
where A — m/m p , ti oss ~ 50Myr for protons or nuclei 
20, 21]. For T jt: r cr i t , particles arc not accelerated to 
Because the classical synchrotron formula used 
in Eq. @ is invalidated when the energy of the emit- 
ted photons exceeds the particle energy [221 ] . we verified 
that j'B'/[mlc 3 /{hZe)} < 1 for T > 1 and t > 10~ 4 s, 
confirming our use of the classical formula in deriving 
Eq. ([5]) for nuclei; electron cooling will be modified by 
QED effects, but we do not consider them here. 

Protons reaching E max lose energy via photomeson 
processes, producing secondary neutrinos. Although the 



maximum neutrino energy is E r 
energy is E mscXtV ~ E max /(20A) 




the average 



; that is, 
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which is maximized for protons. 

Although Eq. ((5} establishes an upper bound on E, 
astrophysical engines may not typically drive outflows 
with t = ti ss ~ 50 Myr. Indeed, unsteady relativistic 
flows develop internal shocks, dissipating their energy on 



a scale r = T 2 cSt [26j]. To evaluate Eq. Q in this limit, 
we scale the variability timescale St to the light-crossing 
time of the largest compact objects, super-massive BHs 
of - 10 9 M (t L ~ 2GM/c 3 
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E max (St) ~ 1 x 10 23 eV 
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r crit w 2020 Z 5 / 7 ^- 4 / 7 (10 4 s/^) 1 / 7 (e B /*) 3/14 /- 1/7 . 
For the 10 M Q BHs that result from the collapse of some 
massive stars, tL ~ 10 -4 s, and i? max ~ 8.6 x 10 21 eV for 
protons. Note the dependence on A and Z; for Z oc A, 
EmAfit) A 6 / 7 , so that is « 30 times higher for 

iron nuclei than for protons. 

UHE neutrinos may also come directly from the de- 
cay of unstable particles in the flow (e.g., 7r ± ). Re- 
quiring t' acc < tj' , where r is the particle decay time, 
t > 8 x 10- 4 s(<5Vl0 4 s) 4 / 7 (e B ^ 8 M 5 / 4 *) 1/7 is required 
for acceleration. Thus, [ir* 1 , rr^, and may be directly 
accelerated, but only for St <C 10 4 s, implying that it 
is impossible to exceed the limit of Eqs. (O and flSJ by 
directly accelerated unstable particles. 

Extreme Beaming. — All of the above expressions are 
valid for the case of a beamed flow with Sr' < r' ± . The 
opposite limit requires that V > 1 / (2\£ /1 / 2 ), and for parti- 
cles to be accelerated one requires R' L < r' ± . These con- 
ditions lead to a unique i? max and a minimum Lorentz 



s max (t,r>i/2* 1 / 2 )^z 
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E 



pi 



Z 10 27 eV 



(9) 



and r min ~ 3 x 10 6 (Z 5 //) 1/2 4 /4 ^4- 2 (W0 1/2 , where 
E p \ is the Planck energy and a is the fine structure con- 
stant. Equation (j9j) is identical to the bound for loss- less 
electromagnetic acceleration Q. We have not empha- 
sized this case here since r m ; n is very high, implying a 
solid angle for the flow of < 10~ 12 sr. 

Discussion. — Our primary result is that for maxi- 
mally efficient acceleration in relativistic shocks, driven 
by maximally luminous astrophysical engines, the maxi- 
mum particle energy is given by Eq. ([5]). For extremely 
highly beamed flows (likely unphysically so), the bound 
is given by Eq. (|9]). These bounds yield maximum and 
average secondary neutrino energies that are within reach 
of current and upcoming experiments. 

Figure [T] shows experimental upper bounds on the 
UHE neutrino intensity above 10 19 eV from ANITA (27|, 
WSRT]28], and FORTE [13], and projected bounds for 
SKA and LOFAR The Waxman-Bahcall flux is 

shown for reference (dashed line) [3(| ■ The vertical solid 
lines show the bounds on energy derived here assum- 
ing A = Z = f = €B = l- Lines I and II show 
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FIG. 1: Neutrino intensity versus energy. Current (ANITA, 
WSRT, FORTE) and projected (LOFAR, SKA) experimental 
upper bounds are indicated. The upper limits on accelerated 
particle energy derived here are indicated by the solid vertical 
lines I-V. Lines III and V (heavy) show the absolute upper 
bounds for spherical and highly-beamed flows (Eqs. [S]and[§]), 
beyond which no particle can be accelerated by a relativistic 
shock launched from any astrophysical engine. Lines II and 
IV mark the average secondary neutrino energies expected for 
these two limiting cases. Line I corresponds to Eq. (|5J. 



E max ,„(5t = 10 4 s) (Eq.© and E maXjV (t = t loss ) (Eq.©. 
The heaviest vertical solid line III marks E max (t — t\ oas ) 
(Eq. [5]), the absolute upper bound on UHE neutrino en- 
ergy for spherical shocks (\P = 1). For acceleration to 
line III, the secondary neutrino spectrum drops sharply 
at E maXtV ~ E max _/20 (line II), as does the heavy dotted 
curve (see below). The absolute bound for beamed flows 
is also shown (line V, heavy dashed; Eq. O, as is E niax _ u 
for such particles (line IV). 

Figure [1] shows that the upper bounds presented here 
are testable. If the primary UHECRs that produce the 
yet- to-be-measured UHE neutrino flux are in fact acceler- 
ated by relativistic shocks, we predict a dramatic down- 
turn in the UHE neutrino flux at or below -E max , v — 
5 x 10 22 eV (line II; Eq. [HI, and most likely well below 
~ 10 22 eV, since our supposition is that actual astro- 
physical engines fail to reach L max , have eb <C 1, and 
may not have Y — r cr i t , or t = t\ oss (Eq. [5]). In turn, the 
observation of a cutoff in the UHE neutrino spectrum 
will put limits on the sources' overall power, since L is 
directly connected with -E max Q. In particular, we can 
invert the expressions above to provide a lower limit on 
the sources' luminosity, given the observation of a single 



UHECR of energy E 21 = E/10 21 eV: 
L^xlO^ergss" 1 ^ 2 * f 



ZA 2 e B V St 



1/2 



(10) 



Importantly, the observation of a single neutrino above 
B max (line III) immediately rules out spherical shock ac- 
celeration as the origin of the primary particle's high en- 
ergy. Indeed, since even for a mono-energetic primary 
proton spectrum at E max (III), the vast majority of the 
secondary neutrinos will have £ mHj „ (II), the detection 
of neutrinos above -E max , u strongly indicates that the pri- 
mary was not accelerated by a relativistic shock. The 
remaining possibilities are that (1) the flow was exceed- 
ingly highly beamed (CI < KT 12 sr), (2) QED effects 
suppressed synchrotron cooling (requiring t < 10~ 4 s at 
L = £ max ), (3) the primary particle was produced by an 
unknown class of near-maximal loss-less electromagnetic 
accelerators, as in [7], or (4) the detected neutrino was 
generated by the decay of a super-massive primary from 
beyond the Standard Model [l[lo|]- 

BH-BH Mergers as UHECR Sources?— The fact that 
L mal (Eq. [3]) is only thought to be approached in gravity 
waves during BH-BH mergers motivates a consideration 
of these events as the primary source of the UHECRs. 

If we assume that the galaxy merger rate at z = (~ 
0.03 Gyr -1 per galaxy; [31]) corresponds to the merger 
rate of Mbh ~ 10 8 M Q BHs, that 7?cr-Mbhc 2 is radiated 
in UHECRs per merger, and that the local number den- 
sity of large galaxies like the Milky Way is 10 -2 Mpc -3 , 
by comparing the volumetric power produced in CRs 
with the inferred UHECR production rate per comoving 
volume (rs 0.7 x 10 44 ergs yr _1 Mpc -3 ; |30() one finds 
that rjcK ~ 10~ 6 is required if such mergers are to dom- 
inate UHECR production. Assuming an E~ 2 UHECR 
spectrum to E max (t = t\ oss ), in Figure 1 we sketch the 
expected neutrino spectrum (dotted curve). Taking the 
local GZK volume to be ~ (100 Mpc) 3 , the local merger 
rate is of order 0.3 Myr -1 , implying that ~ 100 BH-BH 
merger events would contribute to the observed UHE- 
CRs. Although such particles would be expected to trace 
the large-scale distributionof galaxies on the sky, as has 
been claimed by Auger 



32|, the low source density and 



high energy release per event may be ruled out if the 
UHECRs are protons [H (see Ref. [H]). 

Stellar-mass BH-BH mergers might also contribute to 
the UHECR budget, but given the inferred rate per 
galaxy (~ 10 -8 -10 -6 yr _1 0) the UHECR energy pro- 
duction rate per volume is smaller than for supermassive 
BH-BH mergers by ~ 10 3 — 10 5 . Nevertheless, since these 
sources are transient and potentially nearby, they could 
in principle contribute to the observed UHECR budget. 

Can BH-BH mergers accelerate CRs? Simulations 
have argued that the efficiency of production of electro- 
magnetic waves in the vicinity of BH-BH mergers may be 
large, but only for ambient magnetic field strengths far 
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in excess of those thought to accompany accretion [36| . 
Thus, £b in our expressions is the primary unknown in 
forwarding BH-BH mergers as the source of UHECRs. 

Consequences of Maximal Shocks. — If any astrophys- 
ical engine, be it merging BHs or otherwise, reaches 
L m£0i , with r = r cr it, there are observational conse- 
quences beyond the direct detection of UHECRs or neu- 
trinos. In particular, the particles will produce syn- 
chrotron radiation with observed photon energy -E S ynch = 
Th(3/4ir)Y' 2 (ZeB'/mc) ~ 20(10 4 s/St) 1 / 7 TeV. The 
power radiated per particle is ~ E/2irT 2 5t. 

Again considering BH-BH mergers and spherical 
shocks, the total received synchrotron power is L S ynch ~ 

ry C RMBHC 2 /(2^t) ~ 3 x 10 51 (10 4 s/<Si)(WlO- 6 ) ergs 
s" 1 for Af B H = 10 s M . If the particles have an E~ 2 
spectrum extending to E ma-X (5t = 10 4 s), vL s J ndcy ~ 
^ S ynch/(21n(£; max /mc 2 )) - 5 x lO^ergss^ 1 at E synch 
with vL v oc v 1 ! 2 for E < E sync h. The observable rate 
of such events would be ~ 10 yr _1 , and the flux would 
be ~ 5 x 10 -8 GeV s~ 4 cm -2 at a luminosity distance 
of ~ 7 Gpc. Such an event could have the appearance of 
a long-duration gamma-ray burst (GRB), and would be 
observable by the Fermi and Swift satellites. Whether or 
not a sub-class of long-duration GRBs might be ascribed 
to such a mechanism is unknown, but unlikely. Simi- 
lar statements could be made for the analogous events 
from lower mass BH-BH mergers, and the association 
with short-duration GRBs. More likely, the lack of such 
events may limit t]cr to be <C 10 -6 , ee to be <C 1, or 
r / r crit in BH-BH mergers. 

Conclusion. — We have derived robust, experimentally 
testable upper limits to the energies of primary UHECRs 
accelerated by relativistic shocks driven by astrophysical 
engines. For such an engine to accelerate particles to 
£ max ~ 10 24 eV (Eq. it must have L = L max , and 
the relativistic outflow it drives must reach T = r cr it, 
with ee = 1. For highly beamed flows, E max ~ 10 27 eV 
(Eq. [9|), with the additional requirement that the solid 
angle for the flow must be < 10 -12 sr, likely unattainable 
in real astrophysical engines. If any of these conditions 
is not met, E max is not reached. Particles reaching i? max 
suffer photomeson or photodisintegration losses, produc- 
ing secondary neutrinos with average energy ~ 20 times 
lower than E max for protons, observable from across the 
universe. Thus, if a single UHE neutrino is detected 
above -©max it will strongly imply the existence of an 
unknown class of loss-less maximal electromagnetic ac- 
celerators , or physics beyond the Standard Model. 
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